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ABSTRACT 

Due to the increasing heterogeneity and deployment den¬ 
sity of emerging cellular networks, new flexible and scal¬ 
able approaches for their modeling, simulation, analysis and 
optimization are needed. Recently, a new approach has 
been proposed: it is based on the theory of point processes 
and it leverages tools from stochastic geometry for tractable 
system-level modeling, performance evaluation and optimiza¬ 
tion. In this paper, we investigate the accuracy of this 
emerging abstraction for modeling cellular networks, by ex¬ 
plicitly taking realistic base station locations, building foot¬ 
prints, spatial blockages and antenna radiation patterns into 
account. More specifically, the base station locations and 
the building footprints are taken from two publicly available 
databases from the United Kingdom. Our study confirms 
that the abstraction model based on stochastic geometry is 
capable of accurately modeling the communication perfor¬ 
mance of cellular networks in dense urban environments. 

Categories and Subject Descriptors 

C.2.1 [Network Architecture and Design]: Wireless com¬ 
munication; C.4 [Performance of Systems]: Modeling 
techniques; G.3 [Probability and Statistics]: Stochastic 
processes; 1.6.4 [Simulation and Modeling]: Model Val¬ 
idation and Analysis; 1.6.5 [Simulation and Modeling]: 
Model Development. 

General Terms 

Networks, Modeling, Simulation, Performance Evaluation. 

Keywords 

Cellular networks, stochastic geometry, point processes. 

1. INTRODUCTION 

Heterogeneous ultra-dense cellular networks constitute an 
enabling architecture for achieving the disruptive capabili¬ 
ties that the fifth generation (5G) of cellular networks is ex¬ 


pected to provide pp. Modeling, simulating, analyzing and 
optimizing such networks is, however, a non-trivial problem. 
This is due to the large number of access points that are ex¬ 
pected to be deployed and their dissimilar characteristics, 
which encompass deployment density, transmit power, ac¬ 
cess technology, etc. Motivated by these considerations, sev¬ 
eral researchers are investigating different options for mod¬ 
eling, simulating, mathematically analyzing and optimizing 
these networks. The general consensus is, in fact, that the 
methods used in the past for modeling cellular networks, 
e.g., the hexagonal grid-based model [2], are not sufficiently 
scalable and flexible for taking the ultra-dense and irregular 
deployments of emerging cellular topologies into account. 

Recently, a new approach for overcoming these limitations 
has been proposed. It is based on the theory of point pro¬ 
cesses (PP) and leverages tools from stochastic geometry 
for system-level modeling, performance evaluation and opti¬ 
mization of cellular networks [3]. In this paper, it is referred 
to as the PP-based model. Unlike its grid-based counterpart, 
the locations of the base stations (BSs) are not assumed to 
be regularly deployed, but they are assumed to be randomly 
distributed according to a PP. This approach, due to its 
mathematical flexibility for modeling heterogeneous ultra- 
dense cellular deployments, has been extensively used in the 
last few years and it is gaining exponential prominence in 
the scientific community. The interested reader is referred 
to A-8, for the latest achievements in this field of research. 

Despite that, the experimental validation of the PP-based 
abstraction for modeling cellular networks has remained elu¬ 
sive to date. This is especially true for modeling macro cel¬ 
lular BSs, whose deployment is, usually, not totally random. 
A few researchers have tried to justify the PP-based model 
by using empirical data for the locations of the BSs, e.g., 
®, [a, m- These studies have confirmed the potential ac¬ 
curacy and the usefulness of the PP-based model. They, 
however, are based on a small set of data and on simplify¬ 
ing modeling assumptions. Notably, they do not account for 
the footprints of the buildings and rely on simplified chan¬ 
nel models, where line-of-sight (LOS) and non-line-of-sight 
(NLOS) propagation conditions that originate from the pres¬ 
ence of buildings are neglected. This is mostly due to the 
inherent difficulties in obtaining accurate data related to the 
locations of the BSs and of the buildings in urban areas m 
The importance of modeling LOS and NLOS propagation 
conditions, however, has recently been emphasized in sev¬ 
eral papers, e.g., m and references therein. 


In this paper, we investigate the accuracy of the PP-based 
abstraction for modeling cellular networks with the aid of ex¬ 
perimental data. We explicitly take realistic BS locations, 
building footprints, LOS/NLOS channel conditions and an¬ 
tenna radiation patterns into account. More specifically: i) 
the locations of the BSs are taken from a large database 
made available by OFCOM, the independent regulator and 
competition authority for the United Kingdom (UK) com¬ 
munications industries US; and ii) the footprints of the 
buildings are taken from a large database made available 
by Ordnance Survey, the Britain’s mapping agency offering 
the most up-to-date and accurate maps of the UK [14]. Our 
extensive study highlights that the PP-based model is ca¬ 
pable of accurately predicting the performance of cellular 
networks in dense urban environments. It also shows, how¬ 
ever, that their performance highly depends on the channel 
models and the antenna radiation patterns being used. 

In addition, this paper provides another important contri¬ 
bution: we introduce tractable approximations for modeling 
empirical LOS/NLOS conditions due to blockages and prac¬ 
tical antenna radiation patterns, which facilitate the system- 
level simulation and the mathematical modeling of cellular 
networks. The relevance of the proposed approximations is 
twofold: i) simulation, analysis and optimization of ultra- 
dense cellular networks turn out to be easier, since neither 
the actual building footprints nor the actual antenna radia¬ 
tion patterns are needed and, more importantly, ii) they offer 
a transparent interface between telecommunication opera¬ 
tors/manufactures and governmental agencies on one side, 
as well as academic and research organizations on the other 
side, since accurate data on the locations of the BSs and 
on the footprints of the buildings, which often constitutes 
sensible and confidential information, does not need to be 
explicitly provided. By using the proposed approximations, 
this data can be provided in a modified format, which allows 
the former organizations not to unveil sensible and confiden¬ 
tial information on actual network and city deployments. In 
this paper, we prove that the proposed approximations are 
accurate enough for modeling cellular networks. In a com¬ 
panion paper, we have recently verified that the proposed 
approximations are suitable for mathematical modeling and 
for computing tractable utility functions for system-level op¬ 
timization as well. Even though a detailed discussion of 
mathematical modeling is outside the scope of this paper, 
the interested reader may find preliminary results in [15] . 

The remainder of this paper is organized as follows. In Sec¬ 
tion [2] the system model is introduced. In Section [3] the 
proposed approximations for facilitating system-level simu¬ 
lations and mathematical analysis are presented. In Section 
[4] the accuracy of the PP-based abstraction model and of 
the proposed approximations is substantiated with the aid of 
experimental data. Also, main findings and takeaway mes¬ 
sages are discussed. Finally, Section 0 concludes this paper. 

2. SYSTEM MODEL 
2.1 Base Stations Modeling 

In order to test the accuracy of the PP-based model for the 
locations of the BSs, we use experimental data from two 
actual deployments of BSs that correspond to the cities of 
London and Manchester in the UK. As mentioned in Section 
[U this data is obtained from OFCOM [12] . More specifically, 
we consider the BSs of two telecommunication operators: 02 


Table 1: BS statistics from OFCOM - The city of 
London ( A = 4 km 2 ). 



02+Vod. 

02 

Vod. 

Number of BSs 

319 

183 

136 

Number of rooftop BSs 

95 

62 

33 

Number of outdoor BSs 

224 

121 

103 

Average cell radius (m) 

63.1771 

83.4122 

96.7577 


Table 2: BS statistics from OFCOM - The city of 
Manchester (A = 1.8km 2 ). 



02+Vod. 

02 

Vod. 

Number of BSs 

37 

16 

22 

Number of rooftop BSs 

25 

12 

14 

Number of outdoor BSs 

12 

4 

8 

Average cell radius (m) 

125.925 

191.492 

163.305 


and Vodafone. The empirical data is summarized in Tables 
[l] and [2] for London and Manchester, respectively. The fol¬ 
lowing terminology and notation are used: i) “rooftop BSs” 
is referred to the BSs that lay inside a geographical region 
(polygon) where a building is located; ii) “outdoor BSs” is 
referred to the BSs that lay in a geographical region where 
no buildings are located. Information about the locations of 
the buildings is provided in Section [2.21 and iii) A denotes 
the area of the geographical region under analysis. From the 
number of BSs (A f) and A, the density of BSs is obtained as 
Abs = M/A. Accordingly, the average cell radius shown in 
the tables is computed as R c = \J 1/(7tAbs) HU- It is worth 
mentioning that 02 and Vodafone share one BS in the con¬ 
sidered region of the city of Manchester. As expected, the 
data in Tables [T] and [2] confirms that a denser deployment 
of BSs is available in London. 

To study the impact of network densification and the poten¬ 
tial gains of sharing the BSs between telecommunication op¬ 
erators, two scenarios are investigated. In the first scenario, 
the BSs of 02 and Vodafone operate at different frequencies, 
thus they do not interfere with each other. This is equiv¬ 
alent to having just one telecommunication operator in the 
region of interest. Hence, only the BSs of one telecommuni¬ 
cation operator are accessible to the typical mobile terminal 
(MT). In the second scenario, 02 and Vodafone share the 
BSs and they operate at the same frequency. So, a denser 
deployment of BSs is available in the region of interest and 
the BSs of both telecommunication operators are accessi¬ 
ble to the typical MT. Furthermore, full frequency reuse for 
the BSs of the same telecommunication operator and a sat¬ 
urated load traffic model are assumed. This implies that, 
with the exception of the serving BS, all the accessible BSs 
at a given frequency act as interferers for the probe MT. 

As far as the PP-based model for the locations of the BSs is 
concerned, we assume that the BSs are distributed according 
to a Poisson PP (PPP). It is known, in fact, that the PPP 
results in a pure random deployment and that it is more 
mathematically tractable than any other PPs available in 
the literature [10]. In other words, it corresponds to a worst 
case scenario for testing the accuracy of the PP-based model. 

2.2 Buildings Modeling 

To take realistic blockages into account, i.e., LOS and NLOS 
propagation conditions due to the locations and the shapes 




















Figure 1: London case study. Dense urban environment, where the 55.9% of the area is occupied by buildings. 
Horizontal and vertical axis provide distances expressed in meters, (a) Entire region under analysis, (b) 
Magnification of a smaller region, (c) Google map view of (b). (d) Satellite view of (b). 



Figure 2: Manchester case study. Dense urban environment, where the 44.46% of the area is occupied 
by buildings. Horizontal and vertical axis provide distances expressed in meters, (a) Entire region under 
analysis, (b) Magnification of a smaller region, (c) Google map view of (b). (d) Satellite view of (b). 


of buildings (see also Section [273]), we use experimental data 
corresponding to the actual deployments of buildings in the 
cities of London and Manchester. As mentioned in Section 
[l] this data is obtained from OS [14]. In particular, the 
same geographical regions as those of the locations of the 
BSs in Section o are considered. To make sure that the 
data obtained from the two independent websites of OF- 
COM and OS can be merged together, we have verified their 
consistency with the aid of Google maps for the same areas. 
Figures [T] and [2] provide a graphical representation of the ge¬ 
ographical areas under analysis, by merging the data from 
OFCOM and OS. As far as the buildings are concerned, 
their elevation is not considered, since this data is not avail¬ 
able in the database. Therefore, the analysis of its impact 
is postponed to a future research study. 

2.3 Blockages Modeling 

The presence of buildings in dense urban environments con¬ 
stitute an inherent source of blockages, which results in LOS 
and NLOS links. Modeling LOS and NLOS propagation 
conditions constitute an important requirement for assess¬ 
ing the physical layer performance of transmission schemes 
within the 3rd generation partnership project (3GPP) [16| . 
As far as the PP-based modeling of cellular networks is con¬ 
cerned, however, this is often overlooked due to the mathe¬ 
matical intractability of modeling links having multiple states. 
Only recently, in fact, some attempts to take LOS and NLOS 
links into account have been made, e.g., [15] and references 


therein. In this paper, three approaches for modeling block¬ 
ages are investigated and compared against each other. 

2.3.1 Empirical-Based Model 

Based on the locations of the BSs and on the locations and 
the shapes of the buildings described in Sections 12.11 and 
12.21 respectively, LOS and NLOS propagation conditions can 
be empirically taken into account. In this paper, LOS and 
NLOS links are identified as follows. Let a generic “outdoor 
BS” and a generic MT in the region of interest. The related 
link is in LOS if no building is intersected by connecting the 
BS and the MT with a straight line. Otherwise, the link is in 
NLOS. Let a generic “rooftop BS”, the BS-to-MT links are 
assumed to be in NLOS. This is a simplifying assumption 
used in other papers as well El, which seems acceptable if 
no information on the elevation of the buildings is available. 

2.3.2 3GPP-Based Model 

The 3GPP provides a statistical model for identifying LOS 
and NLOS links in several scenarios. More specifically, the 
approach is probabilistic and it does not need any empirical 
information about the locations of the BSs and of the build¬ 
ings. On the other hand, a BS-to-MT link is assumed to 
be in LOS or NLOS with some probability that depends on 
the BS-to-MT distance. In this paper, we focus our atten¬ 
tion only on the MTs that are located in outdoor. This is 
because differentiating between outdoor and indoor propa¬ 
gation would require to introduce different channel parame¬ 
ters. This generalization is postponed to future research. As 









































a case study, hence, we consider the LOS/NLOS model for 
outdoor MTs in a urban micro-cell scenario, which provides 
a good representation of the geographical regions considered 
for London and Manchester in Tables mu and Figures [HE] 
In mathematical terms, the probability that a link of length 
r, i.e., the BS-to-MT distance is r, is in LOS is [16] : 

p(LOS) ^ _ m j n ^ 18/r - , 1} ^1 — ^ (1) 

where min{-, •} denotes the minimum function. The proba¬ 
bility that the link is in NLOS is pf NLOS ) (r) = 1 — p( L °s) ( r ). 

2.3.3 1-State Model 

As mentioned in Section E31 LOS and NLOS propagation 
conditions are often neglected and all links are assumed to 
be either in LOS or NLOS. This case study is considered 
in this paper as well, in order to better understand the dif¬ 
ferences between 1-state and 2-state blockage models. For 
example, if all the BS-to-MT links are assumed to be in 
NLOS, p( NLOS ) (r) = 1 and p( LOS ) (r) = 0 regardless of the 
link length r. 


2.5 Antenna Radiation Pattern 

In cellular networks, directional antennas are typically used 
for enhancing the received power of the intended link and, 
simultaneously, for reducing the other-cell interference. In 
this paper, we are interested in understanding their impact 
on the accuracy of the PP-based model for cellular networks. 
Two widespread antenna radiation patterns are considered. 

2. 5.1 Omni-Directional Model 

The omni-directional model constitutes the baseline for com¬ 
paring any other antenna radiation patterns. In particular, 
g (•) can be formulated as follows: 

g(0) = 1 V6»e[-7r,7r) (3) 

2.5.2 3GPP-Based Model 

The 3GPP provides a reference antenna radiation pattern 
for system-level simulations, which is formulated as [18] : 

10 i<A^ 3 dB 7 ) if \Q\<ip 

(min) (4) 

10 to if ip < \8\ < n 


2.4 Channel Modeling 

In addition to modeling blockages, we consider a practical 
channel model. In particular, path-loss, shadowing and fast¬ 
fading are considered. Let a generic BS denoted by BSi and 
a generic outdoor MT denoted by MT*. 

As for the path-loss, we consider the physical-based bounded 

model l ^ = k^ s ' > (max {ro, ri,*})“ , where ri,* de¬ 

notes the BS-to-MT distance, S = LOS or S = NLOS if the 
BS-to-MT link is in LOS or in NLOS, denotes the free- 
space path-loss, a*- 5 * denotes the path-loss exponent, and 
ro is a positive constant that avoids the singularity of the 
path-loss model for r^* —>• 0. 

As for the shadowing, we consider that it is distributed ac¬ 
cording to log-normal random variable with mean and stan¬ 
dard deviation (in dB) equal to and to a^ s \ respec- 
tively. In this paper, shadowing is denoted by X) p. 

As for the fast-fading, we consider that the envelope of the 
links in LOS and in NLOS is distributed according to a 
Nakagami-m (with m > 1) and a Rayleigh random variable 
with mean power fl, respectively. In this paper, the enve- 
lope of the fast-fading is denoted by h\ . Fast-fading and 
shadowing are assumed to be independently distributed. 

Thus, the received power at MT*, can be formulated as: 


Pr = 


PtG^X^ 
(max {ro, ri,fc}) Q 


(S) 


( 2 ) 


where Pt is the transmit power of BSi and Gi,* is the an¬ 
tenna gain of the BSi-to-MT* link. If BSi is the intended 
BS of MT*, we assume Gi,* = 1. If, on the other hand, BSi 
is an interfering BS for MT*, the antenna radiation pat¬ 
terns of BS; and MT* are assumed to be randomly oriented 
with respect to each other and to be uniformly distributed 
in [—7r, 7r). In this case, we have Gi,* = Sbs (^Osmt (8k), 
where <;bs (•)> Smt (•) are the antenna radiation patterns of 
BSs and MTs, respectively, and 9i, 6k are the angle off the 
boresight directions of BSs and MTs, respectively. In this 
paper, we assume g ( 6 ) = gss (6) = gMT (6) < 1 for every 6. 
Further details about g (■) are provied in Section E© 


where 6^ 3dB) denotes the 3 dB beamwidth, g (mln ) denotes 
the minimum gain, and ip = 0^ 3dB - ) ^/g( min )/12 denotes the 
angle that corresponds to the main lobe jl8i . 


2.6 Cell Association Modeling 

The typical (probe) MT is assumed to be served by any ac¬ 
cessible BS that provides the highest average received power 
to it. Thus, path-loss and shadowing are both taken into ac¬ 
count for cell association. Fast-fading, on the other hand, 
is averaged out and neglected. This is, in fact, the typical 
operating condition based on 3GPP specifications 1161 . 


Let MT* be the typical MT and the probe link be identified 


(S) 

by the subscript “0”. Let Gq / be defined as follows: 


C$ = min <j Cf S J = 




K (s) (max{ro,r-i,fc}) 






x. 


( s ) 


(5) 


where max{-,-} is the maximum function, <0^ is the PP 
of the BSs in state S, and 1 /Gq* is the highest average 
received power at MT* from any accessible BS whose BS- 
to-MT* link is in LOS if 5 = LOS or in NLOS if S = NLOS. 


From ©, the serving BS of MT* is that corresponding 
to the inverse average received power defined as Go,* = 
min {G< BOS) , G' N fc LOS) }, since it provides the best link. 


2.7 Problem Statement 

Let BSo be the serving BS of MT*. Based on © and ©, the 
signal-to-interference-plus-noise-ratio (SINR) at the typical 
MT, MT*, can be formulated as follows: 


PtG 0 ,*^ S * o) /^ S * o) 


SINR = 


+ E 


Se{LOS,NLOS} E ig4 , 


/ 

fcl^VBSo l ,k 


■(S,S o) 


( 6 ) 


where X< s / o) = (PtG^/C%)1 (g$ > G 0 ,* (So) ) de¬ 


notes the generic interfering term, 1 (■) denotes the indica¬ 
tor function that originates from the cell association, So € 
{LOS, NLOS} refers to the LOS/NLOS state of the BSo-to- 
MT* link, denotes the noise power, and the antenna gain 








of the BSo-to-MTfc link is Go,k = 1 for both omni-directional 
and 3GPP-based antenna radiation patterns. 

In this paper, the performance metric for quantifying the 
accuracy of the PP-based model is the complementary cu¬ 
mulative distribution function of the SINR, since it provides 
complete information on its distribution. In addition, it cor¬ 
responds to the coverage probability of a typical MT as 
a function of the link reliability threshold. Let T be this 
threshold, it can be formulated as follows: 

P C ov (T) = Pr {SINR > T} (7) 


3. TRACTABLE SIMULATION AND MATH¬ 
EMATICAL MODELING 

As mentioned in Section [I] the PP-based model of cellular 
networks has been shown to be tractable under simplifying 
assumptions for the blockages, the path-loss functions and 
the antenna radiation patterns [12], [15) . HZ). In this sec¬ 
tion, we introduce some approximations for incorporating 
realistic models for the blockages and for the antenna radia¬ 
tion patterns in a PP-based modeling framework. To intro¬ 
duce and justify the proposed methodology, in particular, 
we assume that the BSs are distributed according to a PPP. 
The accuracy of the PP-based model is assessed in Section 
[4] The main focus of this section is, in fact, on proposing 
tractable approximations for modeling LOS and NLOS links 
that originate from realistic deployments of buildings, and 
for modeling practical antenna radiation patterns under the 
assumption that the BSs constitute a PPP. 


3.1 Tractable Modeling of Blockages 

Due to the presence of blockages (i.e., buildings), a generic 
BS-to-MT link may be either in LOS or in NLOS with a 
probability that depends on the BS-to-MT distance, the lo¬ 
cations of the BSs, as well as the locations and the shapes 
of the buildings. This dependence on the distance and on 
the network topology makes the simulation and the analy¬ 
sis of cellular networks based on the PP-based model less 
tractable and more time-consuming. In this section, we 
introduce a piece-wise approximation for taking LOS and 
NLOS probabilities into account. The proposed approxima¬ 
tion is optimized from the point of view of the typical MT 
and, thus, the spatial deployments of the BSs and of the 
buildings are explicitly taken into account. This makes it 
suitable for system-level performance evaluation and opti¬ 
mization. In this paper, it is referred to as the “multi-ball 
approximation” of blockages. 


3.1.1 The Multi-Ball Model 

The idea behind the multi-ball approximation lies in re¬ 
placing the actual LOS/NLOS probability of a typical MT 
with an approximated function, which still depends on the 
BS-to-MT distance r but is piece-wise constant as a func¬ 
tion of r. In particular, we split the BS-to-MT distance in 
N + 1 regions, which correspond to N balls whose center 
is at the MT. Let the radii of the N balls be d\ < cfe < 
,...,< dff- Since r is generic, i.e., r G [0,+oo), we have 
d\ > 0 and djv < +oo. The multi-ball approximation for 
the LOS/NLOS probability can be formulated as follows: 


N +1 


71=1 L J 


( 8 ) 


where S G {LOS, NLOS}, do = 0, d/v+i = +oo, l[ a ,6) (•) is 
the generalized indicator function defined as lr a m (r) = 1 
if r G [a, b) and (r) « 0 if r £ [a,b), and gf^ > 0 
denotes the probability that a link of length r G [a, 6) is in 
state S. Since a link can be either in LOS or in NLOS, the 
following constraint need to be satisfied: 

V a (S) - ... - V o (S) - 1 

2^ q ld 0 ,dl] ~ Ndjv.djv+d “ 

Se{LOS,NLOS} Sg{LOS,NLOS} 

(9) 

From ©> it is apparent that the larger the number of balls 
N is, the better the multi-ball approximation is. Simula¬ 
tion and mathematical complexity, however, increase as N 
increases. In the next sections, we introduce a methodology 
for estimating N, d\ < cfe <,...,< djv, and q) ; , for 

[d n ~i,d„\ 

n = 1,2, ...,N and S G {LOS, NLOS} in order to find a 
good trade-off between complexity and accuracy. 

3.1.2 Path-Loss Intensity Matching 
The proposed approach for computing all the parameters in 
© is based on matching the intensities of the path-losses of 
the actual blockage model of interest and of the approximat¬ 
ing blockage model in (|5|). More specifically, the proposed 
approach leverages the displacement theorem of PPPs [T9]. 

Let Ti t k be the link distance from a generic BS, BSi, to 
the probe MT, MT*,. From Section 12.41 the path-loss is 

l (S) (n, k ) = k ( s) (max{r 0 ,r i ,4)“ (S) for S G {LOS, NLOS}. 
By assuming that the BSs are distributed according to a 
PPP of density Abs, it can be proved that the PP of the 

path-losses (max {ro, ?"i,fc}) a< ' , * G j is still a PPP, 
whose intensity measure can be formulated as follows: 

A pl ([0, a;)) = 2 ttAbs V A^ ([0, x)) (10) 

VL " {LOS,NLOS} ^ VL n J 


r+oo 

Ag> ([0, x)) = J Pr|z < ' S) (r) < xj p < ' S ' > ( r)rdr 

r+oo . , > 

= J Pr (max {ro, r})“ < xj p^- 1 (r) rdr 


( 11 ) 


Equation ©3 is a direct application of the displacement 
theorem of PPPs. Further details are available in !19| . 

The intensity measure in ©3 is uniquely determined by the 
probabilities pi LOS i (•) and p( NLOS ) (•) as a function of the 
link distance r. Let A( j a ] / tua 0 ([.,.)) be the intensity mea¬ 
sure that corresponds to the actual probabilities of LOS 
and NLOS. For example, it can be obtained by using the 
empirical- and the 3GPP-based blockage models in Sections 
12.3.11 and 12.3.21 Let A^j? prox ^ ([•,•)) be the intensity mea¬ 
sure that corresponds to the multi-ball approximation in 
©3. We propose to estimate N, d\ < efo <,...,< djv, 
and gff J , for n m 1, 2,..., N and S G {LOS, NLOS} in 

[d„-l,dn] 

©3 by minimizing the following utility error function: 

11In (A^ ([0, x M )j) - In (Ap approx) ([0, x M )j) ||“ (12) 

where j*j| F denotes the Frobenius norm and xm is chosen 
in order to capture the main body of AG“ tual ) ([. ; •)) j ;. e ., it 
is close to zero for x > xm- The rationale of ©3 originates 







a: 


(LOS.3GPP) 


+ 


([ 0 , x)) = E(x- KLOsr^ 03 ) [(1/2 )(s;/klos) 2/q!LOS 'H (x - klo S 18“ los ) 

(624.064 - 36exp (-(x/klos) 1/ “ los / 36) (l 8 + 0e/klos) 17 “ LO3 ) + 18(:t/k;los) 17q,L03 ) H(x- klo S 18“ los )] 
a (nlos,3Gpp) ([0,®)) = E(x- knlos 18 qnlos ) (1/2) ((a:/ KNL os) 17 “ NL03 - 18) 2 - 786.C 
+36exp (-(1/36)( 3 : /knlos) 17 “ nlos ) (l 8 + (x/knlos) 17 “ nl ° 3 )] 


1.064 


A 


(S,MultiBall) 


E(x — K (S Vo ( 


^ E (a; — K^di ( 


,(S) 

, ( (x/K^y /a(S) -d; 

[d n _i,d n 

]\\ / J 

,(S) 

[d n _i,d n 

j {dn - G*n-l) E (x - E 


( 13 ) 


([ 0 ,*)) = ( 1 / 2 ) 9 ^, [xj k ( s) J ^ 

+ ( 1 / 2 )l [ ( 0 S d 1 ] rf ?^ (* - K (S) d? <S) ) + (l/2)g[£ i0 o] (^(:e/k ( s) ) ^ W (* - K (s) d„ S) ^ 

dn-i^j e(x- K (S) d“ii) E (x - K (s) d^ ( 


(14) 


from the fact that, from the point of view of a typical MT, 
the impact of blockages is almost the same if the intensity 
measures are close to each other. The logarithm in m is 
used to make the approximation more accurate. 

In practice, the optimization problem in in can be solved 
by using, e.g., the function Isqcurvefit, which is built-in in 
Matlab. We have verified that the solution is quite stable 
with respect to the choice of the initial point of the search, 
which, then, can be chosen at random. The computation 
of in, however, requires the closed-form expressions of the 
intensity measures Ap£ tua1 ^ .)) and A (appr° x ) ([.,.)). They 
are provided in what follows for relevant case studies. 

Empirical-Based Model. The intensity measure is com¬ 
puted based on the actual footprints of the buildings, ac¬ 
cording to Section [2.3.11 Since LOS and NLOS probabilities 
are not available in closed-form in this case, ED cannot be 
directly applied and the intensity measure cannot be formu¬ 
lated in closed-form. For ease of description, the procedure 
for computing it is discussed in Section [4] 

3GPP-Based Model. The intensity measure is computed 
by assuming the LOS and NLOS probabilities in Section 
12.3.21 More specifically, by inserting (JTJ in GD and comput¬ 
ing the integrals, the intensity measure in m is obtained, 
where E (•) is the Heaviside function and E ( x ) = 1 — E (x). 
The result in m holds under the assumption ro < 18 me¬ 
ters, which is usually satisfied for typical setups. 


3.2 Tractable Modeling of Radiation Patterns 

The antenna radiation pattern plays an important role for 
system-level performance evaluation. Incorporating prac¬ 
tical antenna radiation patterns in system-level simulations 
and mathematical analysis may, however, not be straightfor¬ 
ward, as it usually entails a loss of tractability. This is the 
reason why the omni-directional radiation pattern in Section 
I2.5.1l is usually adopted for system-level analysis. In this sec¬ 
tion, we introduce a tractable approximation for modeling 
arbitrary antenna radiation patterns, which can be readily 
used for system-level analysis and simulation. 

Let g( actual ) (.) be the actual antenna radiation pattern of 
interest. We propose to approximate it by using a multi¬ 
lobe antenna radiation pattern, where the antenna gain is 
constant in each lobe. In this paper, it is referred to as 
“multi-lobe approximation”. In mathematical terms, the ap¬ 
proximating antenna radiation pattern can be formulated as 

^(approx) ^ = fl (MultiLobe) whe re: 


^(MultiLobe) _ 


? (1) 

+2) 


if \o\ < e m 
if 6 » (1) < \e\ < # (2) 

(16) 


g (K) if < \e\ < 7T 


with K denoting the number of lobes and 0 < 0 ^ < ... < 
e ( K -D 

< 7 r being the angles that correspond to the I\ lobes. 


Multi-Ball Model. Similar to the 3GPP-based model, the 
intensity measure can be formulated in closed-form by in¬ 
serting m in ED and by computing the resulting integrals. 
The final result is available in (ED for ro < d\ , which usually 
holds for typical setups. 


1-State Model. The 1-state model can be viewed as a spe¬ 
cial case of the multi-ball model, where N — 0, do = 0, 
di = oo, as well as 5 ^°^ = 1 and g[ 0 N ^° S) = 0 or = 0 

and = 1. Thus, the intensity measure follows from 

ED, and it can be formulated as ( S = LOS or S = NLOS): 

A+state) (|o, J,)) = jr^BS (x j K (S) ) ^ E (x — ( 

( 15 ) 


Similar to ED, the parameters of the multi-lobe approxi¬ 
mation in (1161) . i.e., K and 0 < 6^ < ... < < 71 -, are 

computed by minimizing the following utility error function: 

|iog 10 (<? (actual) ( 9 )) - iog 10 ( 5 (approx) (0)) || 2 f (17) 

Similar to ED, the logarithm is used for guaranteeing a 
better approximation. Since the antenna radiation pattern is 
usually measured in dB, log 10 (-) instead of ln(-) is used. It is 
apparent from ED that the larger the number of lobes K is, 
the better the approximation is. The numerical complexity 
increases, however, as K increases. In Section [D numerical 
examples confirm that a good accuracy can be obtained even 
by using a small number of lobes. 












Table 3: 3-ball approximation of empirical and 3GPP blockage models obtained as the solution of (1121) . 



d\ (meters) 

d 2 (meters) 

ds (meters) 

.(LOS) 

Wil 

(Los) 
q \d, ,d 2 1 

(Loy) 

y\d2,dz] 

(LOS) 

9[da.°ol 

London 

15.1335 

56.5978 

195.7149 

0.7948 

0.3818 

0.0939 

0 

Manchester 

13.2076 

57.8840 

213.3940 

0.7866 

0.4981 

0.1015 

0.0001 

3GPP 

47.7989 

215.9387 

1874.442 

0.9446 

0.2142 

0.0243 

0.0021 


London Manchester 3GPP 




Figure 3: Comparison of the intensity measures 
of empirical- and 3GPP-based models against their 
multi-ball (N = 3) approximation counterpart. 

4. NUMERICAL RESULTS: EXPERIMEN¬ 
TAL VALIDATION 

In this section, we illustrate several numerical examples in 
order to validate the accuracy of the PP-based abstraction 
for modeling cellular networks, as well as to confirm the 
tightness of the proposed multi-ball and multi-lobe approx¬ 
imations for simplifying the simulation and for enabling the 
mathematical modeling of cellular networks. 

Simulation Setup. Unless otherwise stated, the following 
simulation setup, which is in agreement with the long term 
evolution advanced (LTE-A) standard, is assumed. The 
transmit power of the BSs is Pt = 30 dBm; the noise power 
is u 2 n = —174 + 10 log 10 (Bw) + -F d B, where Bw=20 MHz is 
the transmission bandwidth and TTib = 10 dB is the noise 
figure; klos = knlos = (47T/T) 2 is the free space path- 
loss at a distance of 1 meter from the transmitter, where 
v = c/f c is the signal wavelength, c«3x 10 8 meters/sec is 
the speed of light, and f c = 2.1 GHz is the signal frequency; 
ro = 1 meter; the path-loss exponents of LOS and NLOS 
links are olos = 2.5 and onlos = 3.5; the mean and stan¬ 
dard deviation of the shadowing are /tlos = /^nlos = 0 dB, 
ulos = 5.8 dB, and onlos = 8.7 dB; the fast-fading en¬ 
velope of the LOS links follows a Nakagami-m distribution 
with parameters m = 2 and Ulos = 1; the fast-fading enve¬ 
lope of the NLOS links follows a Rayleigh distribution with 
parameter LInlos = 1. Used notation: “02” means that only 
the BSs from 02 are accessible; “Vodafone” means that only 
the BSs from Vodafone are accessible; and “02+Vodafone” 
means that all BSs from 02 and Vodafone are accessible. 

Multi-Ball Approximation of the Blockages. In Fig. QT] 
we study the accuracy of the proposed multi-ball approxi¬ 
mation for modeling spatial blockages, under the assump¬ 
tion that the BSs are distributed according to a PPP. More 
specifically, Fig. [3] compares the intensity measure of the 
empirical-based and of the 3GPP-based models against their 
multi-ball approximations, which are obtained by solving 
the optimization in m- The solution of the optimization is 
summarized in Table [3] by assuming N = 3, which provides 


a good matching accuracy while still keeping the computa¬ 
tional complexity at a low level. The intensity measure of 
the empirical-based model of the blockages is obtained by 
using the following procedure based on the actual locations 
and shapes of the buildings obtained from the OS database. 

Step 1: The geographical regions illustrated in Figs. [Qand 
[5] for London and Manchester are considered. As dis¬ 
cussed in Section m this data is obtained from OS. 
The BSs are generated according to a PPP of density 
Abs, which is chosen according to the data in Tables [T] 
and [5] Outdoor and rooftop BSs are identified. 

Step 2: In the same areas, the MTs are generated according 
to another PPP of density Amt = IOAbs- This choice 
of Amt guarantees saturated traffic conditions, i.e., all 
the BSs have at least one MT to serve based on the 
cell association in Section 12.61 Among all the MTs, 
one MT among those that do not lay in a building 
(outdoor MTs) is randomly chosen as the typical MT. 
Step 3: Let the probe (typical) MT, its distance r and link 
state (LOS or NLOS) with respect to any accessible 
BSs are computed according to Section (2.3.11 
Step 4: Step 2 and Step 3 are repeated several thousands 
of times in order to get sufficient statistical data. From 
this data, two vectors are obtained: a vector contain¬ 
ing the distances whose links are in LOS and a vector 
containing the distances whose links are in NLOS. 
Step 5: From the vectors computed in Step 4, p( LOS ) (r) 
and p( NLOS l (r) are estimated by using, e.g., the hist 
function of Matlab. To this end, a resolution step equal 
to Ar = 1 meter and Mt = 2000 discrete distances are 
considered. Thus, the LOS and NLOS probabilities are 
available for the set of distances rt for t = 1,2,..., Mt, 
where Ar = rt — rt -i = 1 meter. The correspond¬ 
ing LOS and NLOS probabilities are pTOS) and 
p (NLOS) (r t ) for £ = 1,2,..., M t . 

Step 6: Finally the intensity measure of the path-losses is 
computed by using CEB and the following discrete (em¬ 
pirical) approximation of (1111 (for S = {LOS, NLOS}): 

A^([0,x)) 

Mt (18) 
« Ar Pr | k ^ (max {ro, r t }) a< ) < x| p' S ' > (r t ) r t 

t =l 

The empirical intensity measure in © constitute the “ac¬ 
tual” intensity for the left and middle plots in Fig. [3] The 
right plot is obtained by assuming m as the “actual” in¬ 
tensity. All in all, a good accuracy is obtained. 

Multi-Lobe Approximation of the Antenna Radiation 
Patterns. In Fig. []] we study the accuracy of the pro¬ 
posed multi-lobe approximation for modeling realistic an¬ 
tenna radiation patterns. In particular, we test the multi¬ 
lobe approximation for modeling the 3GPP-based antenna 























Figure 4: 3GPP-based antenna radiation pattern in 
Q and its multi-lobe approximation obtained by 
using (1161) and (1 171) . The multi-lobe approximation 
is obtained for K = 4 and the solution of (1171) is 

g (1) = 0.8341, g (2) = 0.2865, g (3) = 0.0334, g w = 0.005, 
# (1) = 16.152°, d (2) = 32.304°, and d (3) = 48.455°. 

London (02) London (02+Vodafone) 






Figure 5: Coverage probability by assuming an 
omni-directional antenna radiation pattern. Three 
case studies are analyzed. 1) “OFCOM, OS”: the 
BSs are obtained from the OFCOM database and 
the buildings from the OS database. 2) “PPP, OS”: 
the BSs are distributed according to a PPP and the 
buildings are obtained from the OS database. 3) 
“PPP, Multi-Ball”: the BSs are distributed accord¬ 
ing to a PPP and the multi-ball approximation in 
Table l3l (London and Manchester) is used. 

radiation pattern in ©• The figure shows that a good ac¬ 
curacy is obtained even though K = 4. The approximation 
error can be reduced by increasing K. In what follows, we 
show, however, that the 4-lobe approximation in Fig. [4] is 
sufficiently accurate for system-level performance analysis, 
while offering a low computational complexity. 

PP-Based Modeling of the BSs. In Figs. EE we study 
the accuracy of the PP-based abstraction for modeling cel¬ 
lular networks, by either considering or not the multi-ball 
and the multi-lobe approximations for modeling blockages 
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Figure 6: Coverage probability by assuming an 
omni-directional antenna radiation pattern. Three 
case studies are analyzed. 1) “OFCOM, 3GPP”: the 
BSs are obtained from the OFCOM database and 
the LOS/NLOS link states are obtained from ([l]). 
2) “PPP, 3GPP”: the BSs are distributed according 
to a PPP and the LOS/NLOS link states are ob¬ 
tained from (fljl. 3) “PPP, Multi-Ball”: the BSs are 
distributed according to a PPP and the multi-ball 
approximation in Table l3l (3GPPI is used. 
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Figure 7: Coverage probability by assuming the 
3GPP antenna radiation pattern in ©. Three case 
studies are analyzed. 1) “A”: the BSs are obtained 
from the OFCOM database and the buildings from 
the OS database. 2) “B”: the BSs are distributed 
according to a PPP and the buildings are obtained 
from the OS database. 3) “C”: the BSs are dis¬ 
tributed according to a PPP, the multi-ball approxi¬ 
mation in Table l3l (London and Manchester) is used, 
and the 3GPP antenna radiation pattern is replaced 
by its multi-lobe approximation in Fig. © 


and antenna radiation patterns, respectively. As discussed 
in Sections O and El the empirical coverage probability 
is obtained by using the locations of the BSs and the foot- 






































































prints of the buildings from the OFCOM and OS databases, 
respectively. The following procedure for computing the em¬ 
pirical coverage probability is used. 

Step 1: The geographical regions illustrated in Figs. [l]and 
[2] for London and Manchester are considered. As dis¬ 
cussed in Section E3 this data is obtained from OS. 
Two case studies for the locations of the BSs are con¬ 
sidered. 1) The BSs are distributed according to their 
actual locations obtained from OFCOM (Figs, [l] [2] 
Tables [l] 0. 2) The BSs are distributed according to 
a PPP whose density is the same as that of Tables|T|[2] 
In both cases, outdoor and rooftop BSs are identified. 
Step 2: In the same areas, the MTs are generated according 
to another PPP of density Amt = IOAbs- This choice 
of Amt guarantees saturated traffic conditions, i.e., all 
the BSs have at least one MT to serve based on the 
cell association in Section 12.61 Among all the MTs, 
one MT among those that do not lay in a building 
(outdoor MTs) is randomly chosen as the typical MT. 
Step 3: Let the probe (typical) MT, its distance r and link 
state (LOS or NLOS) with respect to any accessible 
BSs are computed according to Section [2.3.11 
Step 4: For each link between the probe MT and the acces¬ 
sible BSs, path-loss, shadowing and fast-fading gains 
are generated according to Section 12.41 
Step 5 : Let the probe MT and the accessible BSs, its serv¬ 
ing BS is identified by using m in Section m 
Step 6: The antenna gain of the probe link is set equal to 
one, while the antenna gains of all the other links are 
generated as described in Sections 12.41 and 12.51 
Step 7: The SINR and its associated coverage probability 
are computed by using m and 0 , respectively. 

Step 8: Steps 1-7 are repeated 10 6 times in order to get 
sufficient statistical data. The final coverage probabil¬ 
ity is computed as the empirical mean of the obtained 
10 6 realizations for each target reliability threshold. 

If the LOS and NLOS probabilities are computed based on 
the 3GPP-based blockage model in Section [2.3.21 the same 
procedure is used. The only difference is that the actual lo¬ 
cations of the buildings are ignored and each link state (LOS 
vs. NLOS) is identified, according to 0, only based on the 
distance between the probe MT and each accessible BS. The 
same comment applies to the 1-state model in Section [2.3.31 
In this case, all the links are assumed to be, a priori, either 
in LOS or NLOS. 

All in all, Figs. 15171 confirm the accuracy of the PP-based 
abstraction model, and the tightness of the proposed multi¬ 
ball and multi-lobe approximations in practical scenarios. 

Achievable Performance: Impact of Blockages and An¬ 
tenna Radiation Patterns. In Fig. 0 we study the impact 
of the blockage model on the coverage probability. This fig¬ 
ure highlights the importance of accurately modeling block¬ 
ages. More specifically, the widespread used 1-state model 
provides different results from the more accurate and re¬ 
alistic LOS/NLOS blockage model, which accounts for the 
locations of buildings. Figure [8] points out that the coverage 
probability may be better than that predicted by using the 
1-state model, since the links in LOS result in good probe 
links while the links in NLOS result in less interference. The 
proposed multi-ball approximation turns out to be a useful 
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Figure 8: Coverage probability: impact of blockages 
(omni-directional antennas). Three case studies are 
analyzed. 1) “OS”: the buildings are obtained from 
the OS database. 2) “1-State (N)”: all links are in 
NLOS. 3) “1-State (L)”: all links are in LOS. 
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Figure 9: Coverage probability: impact of the an¬ 
tenna radiation pattern (the buildings are obtained 
from the OS database, i.e., empirical-based block¬ 
ages are considered). Three case studies are ana¬ 
lyzed. 1) “02”: only the BSs from 02 are accessi¬ 
ble. 2) “Vodafone”: only the BSs from Vodafone are 
accessible. 3) “02-(-Vodafone”: all BSs from 02 and 
Vodafone are accessible. 


tool for taking LOS/NLOS propagation conditions into ac¬ 
count at an affordable complexity, eventually leading to the 
mathematical analysis and optimization of cellular networks. 


In Fig. [2] we study the impact of the antenna radiation 
pattern on the coverage probability. This figure highlights 
the importance of accurately modeling the antenna radia¬ 
tion pattern. In particular, it shows that sharing the in¬ 
frastructure between telecommunication operators may not 
result in any improvements of the coverage probability if 






























omni-directional antennas are used. If, on the other hand, 
practical (i.e., directional) antennas are used, e.g., based on 
3GPPP recommendations (see 00, the coverage probability 
may be better due to the denser deployment of BSs. This is 
due to the fact that directional antennas have the inherent 
capability of reducing the impact of interference, which is 
desirable in interference-limited cellular networks. 

5. CONCLUSION 

With the aid of experimental data for the actual locations 
of BSs and for the actual locations and shapes of buildings, 
we have studied the accuracy of the PP-based abstraction 
for modeling cellular networks. This study has highlighted 
that the PP-based model is sufficiently accurate for modeling 
dense urban environments of major metropolitan areas. We 
have observed that accurate models for the blockages and 
for the antenna radiation patterns are needed for obtaining 
reliable estimates of the coverage probability of cellular net¬ 
works. Finally, we have proposed flexible approximations 
for incorporating realistic blockage and antenna models into 
the PP-based abstraction of cellular networks, and have val¬ 
idated their accuracy in relevant scenarios. Based on these 
findings, the PP-based model seems to be sufficiently accu¬ 
rate and tractable for enabling the mathematical analysis 
and optimization of emerging ultra-dense cellular networks, 
which use advanced wireless access transmission schemes. 
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